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Positive frequency Wightman function and vacuum expectation values of the energy-momentum 
tensor are computed for a massive scalar field with general curvature coupling parameter and satis- 
fying Robin boundary condition on a uniformly accelerated infinite plate. The both regions of the 
right Rindler wedge, (i) on the right (RR region) and (ii) on the left (RL region) of the plate are 
investigated. For the case (ii) the electromagnetic field is considered as well. The mode summation 
method is used with combination of a variant of the generalized Abel-Plana formula. This allows 
to present the expectation values in the form of a sum of the purely Rindler and boundary parts. 
Near the plate surface the vacuum energy-momentum tensor is dominated by the boundary term. 
At large distances from the plate and near the Rindler horizon the main contribution comes from 
the purely Rindler part. In the RL region the vacuum energy density of the electromagnetic field is 
■ negative near the horizon and is positive near the plate. 
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I. INTRODUCTION 



The influence of boundaries on the vacuum state of a quantum field leads to interesting physical consequences. 
Well known example is the Casimir effect (see [^]-^ and references therein), when the modification of the zero-point 
fluctuations spectrum by the presence of boundaries induces vacuum forces acting on these boundaries. Another 
interesting effect is the creation of particles from the vacuum by moving boundaries In this paper we will 
. study the scalar vacuum polarization brought about by the presence of infinite plane boundary moving by uniform 
T-H ' acceleration through the Fulling-Rindler vacuum. This problem for the conformally coupled Dirichlet and Neumann 
_ massless scalar and electromagnetic fields in four dimensional spacetime was considered by Candelas and Deutsch . 

■ Here we will investigate the vacuum expectation values of the energy-momentum tensor for the massive scalar field 
[ with general curvature coupling and satisfying Robin boundary condition on the infinite plane in an arbitrary number 

■ of spacetime dimensions. The dimensional dependence of physical quantities is of considerable interest in the Casimir 
effect and is investigated for various geometries of boundaries and boundary conditions (see , [|j , Q , [|j , ||] , |l^] 

O I and references therein). The Robin boundary conditions are an extension of the ones imposed on perfectly conducting 
' boundaries and may, in some geometries, be useful for depicting the finite penetration of the field into the boundary 
J> with the " skin-depth" parameter related to the Robin coefficient . It is interesting to note that the quantum scalar 
field satisfying Robin condition on the boundary of a cavity violates the Bekenstein's entropy-to-energy bound near 
. , certain points in the space of the parameter defining the boundary condition . This type of conditions also appear 
in considerations of the vacuum effects for a confined charged scalar field in external fields p3[ | and in quantum gravity 
" ■ 10, |T5[|, Ip^. Mixed boundary conditions naturally arise for scalar and fermion bulk fields in the Randall-Sundrum 
model|17|| . In this model the bulk geometry is a slice of anti-de Sitter space and the corresponding Robin coefficient 
is related to the curvature scale of this space. Unlike to Ref. ||], where the region of the right Rindler wedge on the 
right from the barrier (in the following we will denote this region as RR) is considered only, we investigate vacuum 
energy density and stresses in both regions, including the one between the barrier and Rindler horizon (region RL 
in the following consideration). Our calculation is based on the summation formula derived in Appendix |A| on the 
base of the generalized Abel-Plana formula jlSj (see also |l^). This allows to extract explicitly the part due to the 
unbounded Rindler spacetime and to present the boundary part in terms of strongly convergent integrals. 

The paper is organized as follows. In the next section we consider the vacuum expectation values for the energy- 
momentum tensor of the scalar field in the region on the right of the uniformly accelerated infinite plane. By using 
the summation formula derived in Appendix ^ these quantities are presented in the form of a sum containing two 
parts. The first one corresponds to the energy- momentum tensor for the Rindler wedge without boundaries and is 



investigated in Sec. III. The second term presents the contribution brought about by the presence of the boundary and 



* E- mail : saharyan@ser ver . phy sdep . r . am 



1 



is considered in Sec. 
are evaluated in Sec 



rv. The corresponding boundary terms for the region between the barrier and Rindler horizon 
^ for scalar and electromagnetic fields. In Appendix ^ we show that the direct evaluation of 
the energy-momentum tensor in D — 1 gives the same result as the analytic continuation from higher dimensions. 
And finally, in Appendix we prove the identity used in Sec. Ill to see the equivalency of our expressions for the 
expectation values of the energy-momentum tensor in the case of the FuUing-Rindler vacuum without boundaries to 
the results previously investigated in the literature for the conformally and minimally coupled scalars. 



II. VACUUM EXPECTATION VALUES AND THE WIGHTMAN FUNCTION 



Consider a real scalar field ip{x) with curvature coupling parameter ( in D-dimensional space satisfying Robin 
boundary condition 



{A + Bn'Vi)Lp{x) = 



(1) 



on a plane boundary uniformly accelerated normal to itself. Here A and B are constants, is the unit normal to the 
plane, is the covariant derivative operator. The Dirichlet and Neumann boundary conditions are obtained from 
(|l|) as special cases. Of course, all results in the following will depend on the ratio A/B only. However, to keep the 
transition to the Dirichlet and Neumann cases transparent we will use the form (ffl). The field equation has the form 



(V,V' + + = 0, (2) 

where R is the scalar curvature for the background spacetime. The values of the curvature coupling Q = and 
C = = ^^fy- correspond to the minimal and conformal couplings, respectively. In the case of a flat background 
spacetime the corresponding metric energy-momentum tensor (EMT) is given by (see, e.g., Q) 

T,k = (1 - 2C)5.(p5fc^ + (2C - \l2)g^kd^^di^ - 2C(^V,Vfe(p + (1/2 - 20n?g,k'p'. 

By using the field equation it can be seen that expression (^) can be also presented in the form 



(3) 



Tik = diipdkif + 
and the corresponding trace is equal to 



C - 4 ) g,:feV,v' 



By virtue of Eq. 



= D{C - Cc)V,;VV' + "^V'■ 
for the vacuum expectation values (VEV's) of the EMT we have 



(0 |r,fc(a;)| 0) = lim V,VkG+ix,x') + 



C - 4 ) 5.fcV/V' 



CV^Vfe 



(0|^2(x)|0), 



(4) 



(5) 



(6) 



where |0) is the amplitude for the vacuum state and we introduced the positive frequency Wightman function: 

G+{x,x')^{OMxMx')\0). (7) 

In Eq. (^ instead of this function one can choose any other bilinear function of fields such as the Hadamard function, 
Feynman's Green function, etc. The regularized vacuum EMT does not depend on the specific choice. The reason for 
our choice of the Wightman function is that this function also determines the response of the particle detectors in a 
given state of motion. By expanding the field operator over eigenfunctions and using the commutation relations one 
can see that 



G+{X,X') ^^Lpa{x)ipl{x'), 



(8) 

satisfying 



where a denotes a set of quantum numbers, {ipa{x)} is a complete set of solutions to the field equation 
boundary condition (|ll|). 

Below we will assume that the plate is situated in the right Rindler wedge, xi > \t\. It is convenient to introduce 
Rindler coordinates (r, ^,x) related with Minkowski coordinates (i, a;i,x) by 



t = sinhr, xi = coshr, 



(9) 
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where x = {x2, ...jXd) denotes the set of coordinates parallel to the plate. In these coordinates a wordline defined by 
= const describes an observer with constant proper acceleration The Minkowski line element restricted to 
the Rindler wedge is 



(10) 



The corresponding metric is static, admitting the Killing vector field d/dr. The FuUing-Rindler vacuum is the vacuum 
state determined by choosing positive frequency modes to have positive frequency with respect to this Killing vector. 

First we will consider the region on the right from the boundary, ^ > a, RR region, where is the proper 
acceleration of the barrier. In Rindler coordinates boundary condition (|l|) for a uniformly accelerated mirror takes 
the form 



A + B 



(p — 0, ^ — a. 



The problem symmetry allows to write the solutions to field equation (j^) as 
Substituting this into Eq. (H) we see that the function satisfies the equation 



(11) 



(12) 



(13) 



imaginary order. It can be seen that for any two solutions to equation (|13|), ^^''(0 and (0^{S,), the following 



The corresponding linearly independent solutions are the Bessel modified functions /±iw(A^) and Ki^{\£^) with the 
imaginary order. It can be se 
integration formula takes place 



(14) 



As we will see below the normalization integrals for the eigenfunctions will have this form. 

For the region ^ > a a complete set of solutions that are of positive frequency with respect to d/dr and bounded 
as ^ — » 00 is 



if^ix) = CKUK)e'''^^"^\ a = {lo, k). 
From boundary condition (^ij) we find that the possible values for uj have to be roots to the equation 



(15) 



(16) 



where the prime denotes the differentiation with respect to the argument. This equation has infinite number of real 
zeros. We will denote them by a; = a;„ = a'„(fc), n = 1, 2, ... arranged in ascending order: w„ < w„+i. The coefficient 
C in ([Lq) is determined by the normalization condition 



(</'Q,<Pa') 



'fa dr f*a' = ^ac 



(17) 



with respect to the standard Klein-Gordon inner product, and the ^-integration goes over the region {a,oo). Using 
integration formula (|lj) with cj)^^^ = (/j^^-* — Ki^{X^) and Wronskian W{K^{z), I^{z)} = l/z, from ( [l7| ) one finds 



(27r)^w, 



hu„ (Aa) 



(27r)^-i i^KU^a) I 



where for a given function f{z) we use the notation 

J{z)^Af{z) + bzf'{z)^Q, b^B/a. 



(18) 



(19) 



From the first equality in Eq. (|T^) it follows that is real and positive. Now substituting the eigenfunctions ( p^ ) 
into (0) for the Wightman function we obtain 
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G+{x,x') = 



dk 



(2^) 



D-1 



'ik{x— : 



(20) 



For the further evolution of VEV's (20) we can apply to the sum over n summation formula ( A5). As a function F{z) 
in this formula let us choose 



F{z) = K,,i\OKrA>^e)e-''^^-^'^ ■ 



(21) 



Using the asymptotic formulas for the Bessel modified functions it can be easily seen that condition (A3) is satisfied 
if a^e'"^""^ ' < Note that this condition is satisfied in the coincidence limit t — t' for the points in the region 

under consideration, > a. With F{z) from (|2l] ) the contribution corresponding to the integral term on the left 
of formula (|A5|) is the Wightman function for the Fulling- Rindler vacuum |0/j) without boundaries: 



G+ix,x') = {OR\ipix)ip{x')\OR) 



(2^) 



D-1 



ik(x— x') 



divsmh{TTio)e-"^^^-^'^KUXC)K;U>^^'). (22) 



To see this note that for the right Rindler wedge the eigenfunctions to the field equation (|2|) are in form (see, for 
instance, (2|, pi) 



Vsirdv; 



(23) 



Substitu ting these modes into formula (^ we obtain Eq. (22). Hence, taking into account this and applying summation 
formula (A5) to Eq. (|2C| ) we receive 

(24) 



G+{x,x') ^ Gl{x,x') + {^{x)^{x'))^'\ 
where the second term on the right is induced by the existence of the barrier: 



{if{x)ip{x')) 



'WW - 



dk 



TT J (27r)^-i 



ik(x— x') 



dLo^^.KUXOKU^a cosh[c.(r - r')]. 



(25) 



and is finite for ^ > a. All divergences in the coincidence limit are contained in the first term corresponding to the 
Fulling - Rindler vacuum without boundaries. The VEV's for the EMT are obtained substituting Eq. (|2j) into Eq. 
(^. Similar to (24) this VEV's are also in the form of the sum containing purely Rindler and boundary parts. Firstly 
we will concentrate to the first term. 



III. VACUUM EMT IN THE RINDLER WEDGE WITHOUT BOUNDARIES 

The Wightman function for the right Rindler wedge given by formula (|2^) is divergent in the coincidence limit 
x' X. This leads to the divergences in the VEV for the EMT. To extract these divergences we can subtract from 
Eq. ( p^ the corresponding Wightman function for the Minkowski vacuum |0m )• By using mode sum formula (|^) and 
the standard Minkowskian eigenfunctions it easily can be seen that the latter may be presented in the form 

GUx,x') = I ^^e"^(---')Xo [x^e+e-m'cosHr-r')) , (26) 

where + ^'^ — 2^^' cosh(r — t') = {xi — x[)^ ^ {t ^ ^'Y ■ Using the integration formula 

Eq. ( p^ can be presented in the form convenient for subtraction: 

G+ (x,x') = -^J ^^-^e'^k(x-x ) ^ ^^gj^i^ _ _ r')]}K,UXOKrUXa- (28) 
As a result for the difference between the Rindler and Minkowskian Wightman functions one obtains 
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Gj^{x, x) — Gj^f{x, x') — J- 



1 r d^-^k 



7r2 J {2n 



i-D-l 



„ik(x-x') 



(29) 



After integrating over directions of the vector k, in particular, for the difference between VEV's of the field square it 
follows from here 



X / dLJe~^^cos[uJ{T^T')]KU^0K^^{>^O■ 



(^'(^)>sub = (Oi^|^(.T)|0;^) ~ (OmI(^(x)IOm) 

-1 



dk k 



D-2 



duje-^-Kt{\0 



2D-2^(D+3)/2Y [ D^) 

Using formula (||), for the corresponding difference between the VEV's of the EMT we receive 



(30) 



(31) 



where we have introduced the following notations 



dg{z) 



dz 



dg{z) 



dz 



1 



-2C- 



2C 



(32) 



dg{z) 



dz 



i5(^)r 



\giz)\\ 1 = 2,3,. ..,D, 



and indices and 1 correspond to the coordinates r and ^. 

We consider first the massless limit. In this limit interchanging the orders of integration, making use of the 
integration formula p3] 



•,D-2 



loiO^ / dkk^Ktiki) 



D + l 



r 



D + l 



and the relations 



D- 1 



/ dkk'^-^KUkmLm - -^^id-2{o. 

Jo 



'dkk^KZ{kO-^^l0^ 



for the purely Rindler parts ( |30| ) and (31) one finds 

/<.2x(fl) ^ i / 

W /sub 2^7r^/2+ir(i?/2) Jo 

V » /sub 2^7r^/2+ir(i:)/2) Jo 
where we have introduced notations 



doj e 



dcucu^e-'"^ 



D - 1 



2 

D- 1 



lUJ 



.Ciu) = -Df^^'iu;) = 1 + D{D- l)(Cc - O/co', 
f^piuj) = -l/D +{D- 1)2(C, - C)/cc>2, i = 2, 3, . . . , D. 



(33) 

(34) 
(35) 

(36) 
(37) 

(38) 



Using the recurrence formula, the gamma function in Eqs. ( |3q) a nd (37) can be led to r(l + iuj) for odd D and to 
r(l/2 + iLj) for even D. The standard formulas for the latters ||25[| allow us to write expressions ( |3^ and (37) in the 
form 
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2\{R} ^ 
sub 



/rpk\(R) 

/sub 



-D+l 



2D-inD/2Y(D/2) 



.D-2 



+ (-1) 



— TT 



1=1 



D-1-21 



2oj 



+ (-1) 



D-1-21 
2u! 



(39) 
(40) 



where Im — D /2 — 1 for even D > 2 and Im — {D — \) /2 for odd D > 1, and the value for the product over I is equal 
to 1 for D — 1, 2, 3. These relations illustrate the thermal properties of the Minkowski vacuum relative to the Rindler 
observer [|o), [||. As we see from Eq. (^), though thermal with temperature T = (27r^) ^, the spectrum of 
the vacuum EMT, in general, has a non-Planckian form: the density of states factor is not proportional to uj^~^duj. 
It is interesting to note that for even values of the space dimension D this thermal distribution for the scalar field is 
Fermi-Dirac type. 

Note that for the case of a two dimensional spacetime we can not directly put D = 1 into the integral (|4^). Now, 
though the factor D —1 in the second term on the right of expression ( ^8| ) for flf\uj) is equal zero, the corresponding 
w-integral diverges and this term gives finite contribution to the energy density. Keeping D > 1 the w-integral in Eq. 
( pO| ) can be expressed in terms of the Riemann zeta function (^ji(x), and in the limit D ^ 1 one obtains 



/rpO\{R) 

V-'O/sub 



<2 



^-12.1' 



(41) 



where Cfl(O) = -1/2. 

The density of states factor is Planckian for the conformally coupled scalar and for D = 1,2,3: 



/sub — 



-D-l 



2D-It,d/2y[d/2) Jo e^^^ + (-1)^ 



rdiag( 1,--^,...,--^ ) , C = Cc 



D = 1,2,3. 



(42) 



This result for £) = 3 case was obtained by Candelas and Deutsch Expression (^2|) corresponds to the absence 
from the vacuum of thermal distribution. This means that the Minkowski vacuum corresponds to a thermal state 
with respect to uniformly accelerated observers |2^ . 

We now treat the case of a massive scalar field. First of all let us present the difference (|2^) between Rindler and 
Minkowski Wightman functions in another alternative form. For this we will follow the procedure already used in 
Ref. P3| for the Green function. Substituting |2q] 



+00 



where 7^ = + C'^ + cosh y, into Eq. (p9[) and integrating over uj we receive 



^j:k(x-x') 



+ 00 



(2^)^ 

Integrating over directions of the vector k with the help of formula 



dy ■ 



(y + T-T'y 



J d^'-'k^'^^^-^'^Fik) = (27r)(^-i)/2^ 



(fc|x-x'|)(^-3)/2 

with (x) being the Bessel function, and using the value for the integral over k , 



dk k''+^Jf,{bk)Ko{jVk^ +m^) = &''TO''+^i^p+i(TOV^^M^), 



one obtains 



where — 



G^{x, x) — G\,j{x, — 



+ 00 



dy 



y2 \27r71 



{D-l)/2 



K(D-i)/2i'm"fi), 



(43) 



(44) 



(45) 



(46) 



(47) 



■'2 + 2^^' cosh(2/ — T + t') + |x — x'l . The similar relation for the Feynman's Green function is 
considered in Ref. |23| . Noting that 71 is just the geodesic separation of the points x — (r, x) and x" = {r'+y, — <^', x) 
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and taking into account the standard expression for the Minkowskian Wightman function, relation (47) can be written 
as 



G+ix,x')=G+ix,x') 



+00 



dy 



y 



Gl,{x,x"). 



(48) 



FoUowing Ref. |2^ ], the second term on the right of this formula we can interpret as coming from an image charge 
density — l/(7r^ distributed on a line (— ^',x) in the left Rindler wedge. 

In the coincidence limit Eq. (p7|) yields 



2™ 



D-l 



(27r)(-D+i)/2 



2+y2 z(D-l)/2 ' 



2m^ cosh(?//2). 



(49) 



This quantity is monotone increasing negative function on ^. Substituting Eqs. and (49) into formula (||) for the 
difference between the corresponding VEV's of the EMT we obtain (no summation over i) 



(To' 



0\(fl) 
sub 



2m^+i 



(27r)(^+i)/2 7o y^ + TT 



dy 



[D-l)/2 



(^) 



1 /sub 



lrpi\{R) 

N-' j / sub 



dy 



2m^+i 



(27r)(^+i)/2 7o y-^+TT-^ ziD+i)/ 



z(D-l)/2 

(D+l)/2{z) , 
,(D+l)/2 V 



D- 



K 



(D+l 



)/2{z) 



2m^+i r°° dy 
(27r)(^+i)/2 Jo 

K{D-l)/2{z) 



+ 



AD-l)/2 



[D-l) 



z{D+l)/2 

)/2W 



(D+l 



ziD + l)/2 



z{D+l)/2 

4Ccosh2(y/2)), 
(4C-l)cosh2(y/2) 



(l-4Ccosh2(j;/2)), 



(50) 
(51) 
(52) 



2,3,, 



Note that near the horizon, in the limit <^ 1, the leading terms of the asymptotic expansion for the vacuum 
EMT components (|50|)-(|5^) do not depend on the mass and coincide with Eq. (^). Hence, near the horizon these 
components diverge as ^ . The VEV's of the EMT for the real massive scalar fields with minimal and conformal 
couplings on background of the Rindler spacetime were considered in Ref. by using the covariant functional 
Schrodinger formalism, developed in Ref. In this paper for the EMT the form (||) is used with C = and C = Cc 
and the corresponding formulas [formulas (2.39) and (2.45) in Ref. ] are more complicated. The equivalency of 
these formulas to our expressions (^o|)-(^^ with the special values ^ = 0, can be seen using the identity 



m2^2 



dy 



7r2 ~ 3y2 K^D_iy2{z) 
(7r2+j;2)3 ;j(D-l)/2 



dy 



(1 — cosh J/)' 



(D-l 



)/2{z) 



[£> + 1- (D-l) cosh y] 



z{D-l)/2 

K(D+i)/2{z) 
z(D+l)/2 



(53) 



where z is defined as in Eg. (^9|). The proof of this identity we give in Appendix |^. 

From expressions (^9|)-(^2|) it follows that for the massive scalar field the energy spectrum is not strictly thermal 



and the quantities (Oa/ \v^\Qm) - {Qn, \v^\Qr) = - {^"^ and (Oa/ Oa^) - (0^ \T^\ Qr) = -{T^{x))[^l do not 

coincide with the corresponding ones for the Minkowski thermal bath, and {T^Y^^, with temperature (27r^)~^ 

(for this statement in a variety of situations see Refs. j29j- and references therein). To illustrate this for the case 
D = 3 we present in Fig. ^the quantities — ('/2^)sub (curve a) and 



(cj2 -(_ 2mcj) 



1/2 

2mLjj) dio 



g(m+w)/T _ I 



T ■ 



1 

2< 



(54) 



(curve b) in dependence of m^. 
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FIG. 1. The D — 3 expectation values for the field square in the case of the Rindler vacuum without boundaries, 
(curve a) and for the Minkowski thermal bath, (,'^{'fi'^)'M (curve b), given by formula (|54|), versus m^. 

The integrals in formulas (^)-(|5^) are strongly convergent and useful for numerical evaluation. In Fig. ^thc results 
of the corresponding numerical evaluation for the vacuum EMT components are presented in the Rindler case without 
boundaries for minimally and conformally coupled scalars. 



0.0005 b 





0.5 



1.5 2 2.5 



FIG. 2. The expectation values for the diagonal components of the energy-momentum tensor, ^^^^(T/)'^^ (no summation 
over i), versus in the case of the Fulling- Rindler vacuum without boundaries for minimally (left) and conformally (right) 
coupled scalar fields in D = 3. The curves a, b, c correspond to the values i = 0, 1, 2 respectively. 

In the limit 3> 1 using the asymptotic formula for the Macdonald function and taking into account that the 
main contribution into the integrals in Eqs. ([50|), ([5l|), and ( |5^ ) comes from the region with small y, in the leading 
order one finds 



lrpO\{R) /rpi\(R) 

V^O/sub 



sub 



^' 2^7r(^+3)/2^(^+i)/2 



' sub 



2, D, — > oo. 



In the massless limit, m — 0, from Eq. (49) for the field product one obtains 



The corresponding formulas for the VEV's of the EMT are written as 

,^oU«)_ n/Tl^(«)- ^^(^) r dy 4Ccosh2(y/2)-l 

\-'0/sub — ''^X-'l /sub 



D+l 



y^ + TT^ cosh^+\y/2) 



r (D+l 



dy l-(l-4C)(I?-l)cosh2(y/2) 



13+1 



/o y + 



cosh^+i(2//2) 



= 2, 



.,D. 



(55) 

(56) 

(57) 
(58) 
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Note that the relation ( |57| ) between the energy density and the vacuum effective pressure also follows from the 
continuity equation for the EMT (see Eq. ( |63|) below). Formulas (|5^)-(|5^) are alternative representations for the 
VEV's ®, (H). 



IV. VEV'S IN THE RR REGION 

In Sec. H we have seen that the Wightman function in the case of the presence of a barrier can be presented in the 
form (|^), where the first summand on the right corresponds to the Fulling- Rindler vacuum without boundaries and 
the second one is induced by the barrier and is given by expression (p5|). In the coincidence limit from this formula 
for the boundary part of the field square we have 

This quantity is monotone increasing negative function on ^ for Dirichlet scalar and monotone decreasing positive 
function for Neumann scalar. Substituting the Wightmann function (24) into Eq. (^ and taking into account Eqs. 
(H) and (|5|) for the VEV's of the EMT in the region ^ > a one finds 

(0|2;'=(:r)|0) = {Qn\Tl'{x)\QR) + {Ttix))^'\ (60) 

where the first term on the right is the VEV for the Fulling-Rindler vacuum without boundaries and is investigated 
in previous section. The second term is the contribution brought by the presence of the barrier at ^ = a: 

where the functions F'^^\g[z)\ , z = 0, 1, 2 are obtained from the functions f'^'^\g{z)] in Eqs. (^) replacing uj iuj: 

F(%{z)] = f(%{z),oj^ioj]. (62) 

It can be easily checked that both summands on the right of Eq. ( |60| ) satisfy the continuity equation T^.^ — 0, which 
for the geometry under consideration takes the form 

^(^ri)-ro" = o. (63) 

By using the inequalities 

liiz) < v/TTT^/^(z), ~K{z) > ^l + oj^/z^K^{z) (64) 

it can be easily seen that i^W[ii:^(z)] > 0, F^^^K^iz)] < for C < and i^W[i\:^(z)] > 0, i = 2, 3, ...,£> for 
C ^ aIjT-i) • cases of Dirichlet and Neumann scalars for the boundary parts of the EMT components these 

leads 

((ro")(^)) <0, >0 (65) 

V / Dirichlet V / Dirichlet 

((ro")(^)) <0, ((Ti)W) >0 (66) 

V / Neumann V / Neumann 

for C ^ ^iid (no summation over i) 

({Tl^'A <0, ({T^y'A >0, z = 2,3,...,D (67) 

V / Dirichlet \ / Neumann 

forC<4^. 

For D — 1 case the corresponding boundary part is obtained from Eq. ( |6l| ) 

(^Tl^)i^)=Jlm' Tdcj^^F^nKUmO], ^ = 0,l. (68) 
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Taking the limit m ^ for the massless field we receive from here 



with 



F^"'^^2lj + 1, F^''^ = -1. (70) 



In this case our assmnption for the absence of imaginary eigenmodes corresponds to the restriction AB < on the 
Robin coefficients. We have obtained result ( ]69| ) by using the analytical continuation from the values D > 1. If we 
want to realize a direct evaluation in D = 1 the prescription described in Sec. |l| has to be modified for the massless 
case. In appendix B we show that the corresponding evaluation leads to the same result. 

Boundary part (^) is finite for all values ^ > a and diverges at the plate surface ^ = a. These surface divergences 
are well known in quantum field theory with boundaries and are investigated for various type of boundary conditions 
and geometries. To extract the leading part of the boundary divergence note that near the boundary the main 
contribution into the w-integral comes from large values of uj and we can use the uniform asymptotic expansions for 
the modified Bcssel functions psf . Introducing a new integration variable k —^ ujk and replacing the Bessel modified 
functions by their uniform asymptotic expansions in the leading order one obtains 

(To°(.))'" - in(.r> - ^g^^^«,;,;>'f ;^«. r (^) , (71, 

''-^f^'S-«„ rr£±lV ,72, 



These terms do not depend on mass and Robin coefficients, and have opposite signs for the Dirichlet and Neumann 
boundary conditions. Note that the leading terms for {Tq{x))^''^ and (^2(2;))^^'' are the same as for the plate in 
Minkowski spacetime (see, for instance, [||). This statement is also valid for the spherical and cylindrical boundaries 
on background of the Minkowski spacetime ||], ]To| . 

Now let us consider the asymptotic behavior of the boundary part (^) for large ^ a. Introducing in Eq. 
( |6l| ) a new integration variable y = XS, and using the asymptotic formulas for the Bessel modified functions for small 
values of the argument, we see that the subintegrand is proportional to (ya/2^)'^". It follows from here that the main 
contribution into the w-integral comes from the small values of uj. Expanding with respect to uj in the leading order 
we obtain 

^'^^^ 2^3^(^-i)/2(l + ,5^o)r(V)ln2(2^/a)' ^ ^ 

where 

A^\mO^ dyy\v^-m^ef''-^^'^F^^\K^{y)%=^. (74) 

If, in addition, one has ^ 1 the integral in this formula can be evaluated replacing the Macdonald function by its 
asymptotic for large values of the argument. In the leading order this yields 

4°)^4')^-2mC4^)^7r(l/4-C)r(^)(mO(^+^)/^e-2"«, » 1, z = 2, 3, • • • , i?. (75) 

For the massless case the integral in Eq. (Q) may be evaluated using formulas (|3^)-(|35|) and one obtains 

Ao - - {D-\YT{D)^ \ 2 j' ^^^> 

For the conformally coupled scalar the leading terms vanish and the VEV's are proportional to ^^^^^ ln^'^(2^/a) 
for ^/a 1. As it follows from Eqs. (^), (p5|), and ( [ts]) in both cases ( [75|) and ( [76|) and in the region far from 
the boundary the VEV's (^^ are dominated by the purely Fulling- Rindler part (the first summand on the right): 

{T^)'^'''^ / {T^)[^l ~ ln"^(2^/a). From Eq. (^ we see that for a given ^ the boundary part tends to zero as a ^ 
(the barrier coincides with the right Rindler horizon) and the corresponding VEV's of the EMT are the same as for 
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the Fulling- Rindler vacuum without boundaries. Hence, the barrier located at the Rindler horizon does not alter the 
vacuum EMT. 

And finally, we turn to the asymptotic a, ^ ^ oo, ^ — a = const. In this limit ^/a 1 and w-integrals in Eq. 
( |6l| ) are dominated by the large uj. Replacing the modified Bessel functions by their uniform asymptotic expansions, 
keeping the leading terms only and introducing a new integration variable v ^ cu/a for the energy density one obtains 

(r")W T^^TTTT / dv- 

X 



A + BV^-. 4C.^ + A^(1-4C)]. (77) 



Converting to polar coordinates, v — rcos9, k = rsinO, we can integrate over angles with the help of the formula 



Introducing the new integration variable y — y/r^ + rn^ we receive 



A + By 
A- By 



4(C - Cc) + 



m2(4C- 1) 



a,^— >oo, ^ — a = const. (79) 



Here we have taken into account that in the limit ^ ^ oo the purely Rindler part (O/jITq jO/j) vanishes and in Eq. (|6£ 
the boundary term remains only. In the same order the evaluation of the components, corresponding to the vacuum 
effective pressures in the plane parallel to the barrier, leads (no summation over i) 

{0\t:\0)-{0\T^\0), t^2,3,...,D, (80) 

The leading term for the component (Ojrj'^lO) vanishes, 

(OlTi^lO) = 0(l/a), a^oo. (81) 



The term to the next order can be found using continuity equation (|63|) and expression ([79|). The expectation values 
([tqI), (^0|), and (81) coincide with the corresponding VEV's induced by a single plate in the Minkowski spacetime 



with the Robin boundary conditions on it, and are investigated in Ref. for the massless field and in Ref. |^ for 
the massive scalar field. 

V. VACUUM STRESSES IN THE RL REGION 
A. Scalar field 

In previous section we have considered VEV's for the EMT in the region ^ > a. Now we will turn to the case of 
the scalar field in region between the barrier and Rindler horizon (RL region), < ^ < a, and satisfying boundary 
condition (^ij) on the boundary ^ = a. For this region we have to take both solutions to equation (^3|) and the 
eigenfunctions have the form (^2|), where now 

</-(€) = Z^X^, Xa) EE if,^(AO - ^|^^/,;.(AC), (82) 



and the coefficient is obtained from boundary condition (p, Here the quantities with overbars are defined in 
accordance with Eq. (p^). As we see, unlike to the previous case, the spectrum for lo is continous. The normalization 
coefficient is determined from condition (|l7|), where integration over ^ goes between the limits (0, a). Using integration 



formula (14) and the representation for the delta function 



S{x) = lim , (83) 

0-— ►oo TTX 
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it can be seen that 



Substituting eigenfunctions (|l^ with the function (^) into the mode sum for the Wightman function one finds 



G+{x,x') = -^ I e"-(— / dusiiM^uj)e-''^^^-^">Z,U^^,Xa)Zl{XC,Xa). (85) 



1 /• d^-^k 
(2^)^ 

The part induced by the barrier can be obtained subtracting from here the Wightman function for the right Rindler 
wedge without boundaries, given by expression (^). To transform the corresponding difference note that 

Zt smh TTUJ ^ li^^ ( Xa ) 

where we have used the standard relation 

= (87) 
This allows to present the boundary part in the form 

Assuming that the function Ii^{Xa) (J_i(j(Aa)) has no zeros for — 7r/2 < argw < (0 < argw < tt/2) we can rotate 
the integration contour over to by angle — 7r/2 for the term with a = I and by angle tt/2 for the term with cr = — 1. 
The integrals taken around the arcs of large radius tend to zero under the condition |^^'| < a^e''^"'^ ' (note that, in 
particular, this is the case in the coincidence limit for the region under consideration). As a result for the difference 
of the Wightman functions one obtains 



(^{xMx'))^") ^ — / ^^e^k(x-x') / d^ll^j^^xOU^acosh[cj{T~T')]. (89) 



'^J (27r)^-i' io U\a) 

As we see this expression differs from the corresponding formula ( p5| ) for the region a < ^ < c« by replacements 
— > K^, lui- In the coincidence limit for the VEV of the field square we receive 

(^'(-)>^'' = „n_, .n'.L.^-1^ rdkk^-' rdu^S^lliXO. < ^ < a. (90) 



This quantity is monotone decreasing negative function on ^ for Dirichlet scalar and monotone increasing positive 
function for Neumann scalar. From Eq. (|6|) the VEV's of the EMT are obtained in the form (^0|) with the boundary 
part 



(91) 



where, as in Eq. (pll ), the functions F*^*^ 
Comparing Eqs. (|g), (|9l|), (H), and (|l 



,g{z)] are obtained from expressions ( |32| ) replacing lo ^ ilu (see Eq. (|62|)). 
we see that the boundary parts for the region < ^ < a can be obtained 
from the corresponding ones for the region ^ > a by the replacements I^^ K^, K^^ 1^^. Note that here the 
situation is the same as for the interior and exterior regions in the case of cylindrical and spherical surfaces on 
background of the Minkowski spacetime (see, Q, (lOj). The reason for this analogy is that the uniformly accelerated 
observers produce world lines in Minkowski spacetime which in Euclidean space would correspond to circles. This 
correspondence extends also to the field equation and its solutions, and modes ( p^ ) are the Minkowski spacetime 
analogue corresponding to cylinder harmonics. 

For D = I case the integral over k in Eq. (|8^) is absent and X = m. The corresponding formulas for the field 
square and vacuum EMT are obtained from Eqs. ( |90| ) and (^ replacing 
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OO 1 

dkk^-^^-, X^m. (92) 



In particular, for the massless D ~ 1 case we have 

(Tt)'"' -^Tdu. (93) 
I"? Jo ^ + -Dw/a 

where 

Fo^°) 2cj - 1, = 1. (94) 

The expressions (^ij) are divergent on the plate surface, ^ = a. The leading terms for the corresponding asymptotic 
expansions are derived by the way similar to that for Eqs. ( |7l| ) and (|7^), and can be obtained from these formulas 
replacing ^ — a by a — ^ and changing the sign for {Tl)^^\ 

Now let us consider the behavior of the boundary part (^l|) for ^ <^ a. Introducing instead of fc a new integration 
variable y = Xa and using the formula for Ii^{y^/a) in the case of small values of the argument we see that the w - 
subintegrand is proportional to e"'^'^'"'^"/^^). It follows from here that the main contribution comes from the small 
values of uj. Expanding over u) we obtain the standard integrals of the form aj"e~^"'"('^/^^)(itJ. Performing these 
integrals one obtains in the leading order 



^ °^ ~ ^^1^ ^ 2^-i7r(^+i)/2r(^)a^-ie2 1n2(2a/0 ^ ' 



(4C-l)i?o(ma)(l + 0(ln-\2a/O)) 



where 



BQ{ma) 



' dyy{v^~m^a^r'^y^^. (97) 



Note that the function BQ{ma) is positive for Dirichlet boundary condition and is negative for Neumann one. For 
D — 1 case one has 



^ °' ^ ^' 2<2in2(2/mO loima) ' 



As we see the boundary part is divergent at Rindler horizon. Recall that near the horizon the purely Rindler part 
behaves as and, therefore, in this limit the total vacuum EMT is dominated by this part. 



B. Electromagnetic field 

We now turn to the case of the electromagnetic field in the region < ^ < a for the case D — 3. We will assume 
that the mirror is a perfect conductor with the standard boundary conditions of vanishing of the normal component 
of the magnetic field and the tangential components of the electric field, evaluated at the local inertial frame in which 
the conductor is instantaneously at rest. As it has been shown in Ref. |^ the corresponding eigenfunctions for the 
vector potential Af^{x) may be resolved into transverse electric (TE) and transverse magnetic (TM) (with respect to 
^-direction) modes: 

A(<')(^^-^ iO,0,-ik3,ik2)f^°Hx), for a = , . 

^ ' \ (-ec)/aC,zc^/^, 0,0)^(1) (:e), for a-1' ^ ' 

where k — (^2,^3), cr = and cr = 1 correspond to the TE and TM waves respectively. From the perfect conductor 
boundary conditions on the vector potential we obtain the corresponding boundary conditions for the scalar fields 

(^(")(^): 

^(0)(^)^0, ^^^ = 0, e-a. (100) 
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As a result the TE/TM modes correspond to the Dirichlet /Neumann scalars. In the corresponding expressions for 
the eigenfunctions vlL'^^ (x) the normahzation coefficient is determined from the orthonormahty condition 



OS w 



From here for the corresponding Diriclet and Neumann modes we have 



V sinhc 



zkx— zcjr 



(101) 



(102) 



where k''^{z) = Ki^{z), Kl^\z) — dKii^{z)/dz and the same notations for the function Iiuj{z). The VEV's for the 
EMT can be obtained substituting eigenfunctions (BOf) into the mode sum formula 



'(1 



(103) 



with the standard bihnear form for the electromagnetic field EMT: 



1 



Tt{A^{x),A^{x)} = — -FuF'^ + -S^F,^F 



1 



where Fik = dA^/dx^ — dAi/dx^ is the field tensor. Substituting eigenfunctions 
finds 



{G\T^{x)\0) 



4-^-3 



E 

CT = 0,1 



dk 



dui sinh(aj7r)/4*^ 



(104) 

) into the mode sum ( |103D one 



lt\ka) 



■huikO 



(105) 



Here for a given function g{z) the following notations are introduced 



/i2 [5(^)1 



dg{z) 



dz 
dg{z 



dz 



1 + i5(^)r, 



i.9wr 



(106) 



f(3) 



It can be easily checked that components (105) satisfy zero trace condition and continuity equation (|63|). The way 
for subtraction from these quantities the parts due to the Fulling- Rindlcr vacuum without boundaries is the same as 
for the scalar field, given above in this section. This yields to the following result 



47r2 



dk t' 



duj 



K^ka) , Kiika) 



I^{ka) lUka) 



F^'2[UkO], 



where 



dLu(uj^ +U}) ^ ( 1 1 1 
diag(l,--,--,-- 



(107) 



(108) 



are the VEV's for the FuUing-Rindler vacuum without boundaries and the functions Fctk [<7(z)] are obtained from 
Eq. (|106|) replacing uj icu: Fcm [g{z)] — /im [giz), lu ~* iuj]. Note that the w-integral in Eq. (108) is equal to 11/240. 



Comparing the boundary part in Eq. ( |107D to the corresponding expression for the RR region, £, > a, derived in 
we see that, as in the scalar case, these results arc obtained from each other by the replacements I^^ K^, I^. 
By using inequalities (^J) we conclude that 



K^{z) , K{z) 



IM) 



I'M) 



< 0. 



(109) 
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This, together with the observation that Iuj{z) is monotone increasing function, immediately shows that the com- 
ponents (T!^) z = 2,3 are negative and monotone decreasing on ^. Further from the first inequahty (|6^ ) one has 

Fem[Iu:ikS,)] > and it can be seen that this function is monotone increasing on ^. Hence, the component (T^)''''') is 



we see that Fim^ 



Iuj{k£,)] > 0, and as a result one 



positive and monotone increasing. And, at last, by using Eq. (|f 

obtains that the boundary part of the energy density is positive, (j'0) (^) > Q. 

Now let us consider the asymptotics of the boundary part in Eq. ( |l07| ) near the barrier, ^ — s- a. In this limit the 
corresponding expressions are divergent. It follows from here that the main contributions into the VEV's are due to 
the large u. Rescaling the integration variable k — *■ uk and replacing the Bessel modified functions by their uniform 
asymptotic expansions in the leading order one finds 



307r2o(a-^)3' 



i = 2,3, 



(110) 



l\(b) 



607r2a2(a-C)2' 



^ -> a - 0. 



Note that the corresponding asymptotics for the region ^ > a are given by the same formulas Q . 

Now we turn to the asymptotics in the limit ^ —^ 0. Replacing the function /^^(fci^) by its asymptotic, and noting 
that the main contribution into the w-integrals come from the region with small lu, in the leading order we have (no 
summation over i) 



87r2a4 ln(2a/^) Jq 



dyy^ 



h{y) h{y) 



(111) 



In this formula the value for the integral is equal to -1.326. As we see, unlike to the case for the scalar field, the 
boundary part of the electromagnetic vacuum EMT tends to zero as ^ ^ 0. 

And finally, let u s con sider the leading terms in the vacuum EMT in the limit a, ^ oo, ^ — a ~ constant. 
Introducing in Eq. ( |107| ) a new integration variable ka we see that in this l imit the corresponding terms have the 
same form as in the case ^ — s- a, a = constant and are given by expressions ( 110 ). Note that in this limit the total 
vacuum EMT is dominated by the boundary part. 

In Fig. ^ we have plotted the electromagnetic vacuum energy density, a*{TQ) (left), and perpendicular pressure, 
—a^{Tl) (right), in dependence of ^/a, generated by a single perfectly conducting plate uniformly accelerated normal 
to itself. For the plate ^ = a. The dashed curves present the corresponding regularized quantities for the FuUing- 
Rindler vacuum without boundaries, given by the second summand on the right of formula (lOJ). The full curves are 
for the boundary parts, second term on the right of formula (107). In addition to the region < ^ < a, investigated 
in this section, we have also included the region a < ^ < oo. The formulas for the boundary part in the latter case are 
derived in Ref. js) and, as it have been mentioned above, are obtained from the ones considered here by replacements 
Ki^, Ii^. From Fig. ^we see that the vacuum perpendicular pressure is always negative. 
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FIG. 3. The expectation values for the boundary part of the energy density, a'*(ro) 
tive pressure, — a*(r]^)'''' (right), multiplied by a*, for the electromagnetic field versus ^/a 
corresponding regularized quantities for the Fulling- Rindler vacuum without boundaries [see Eq. (108)] 



(left), and perpendicular effec- 
Th e da shed curves present the 
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The total energy density (as a sum of purely Rindler and boundary parts) is always negative in the region a < ^ < oo. 
For the region < ^ < cxd the total energy density is negative near the horizon and is positive near the plate. As a 
result for some intermediate value of ^ the energy density vanishes. 



VI. CONCLUSIONS 



Quantum field theory in accelerated systems contains many of the special features produced by a gravitational field 
without the complications due to a curved background spacetime. In this paper we have considered the VEV's of the 
EMT induced by uniformly accelerated plane boundary in the FuUing-Rindler vacuum. This problem for the region 
on the right of the plane has been previously studied by Candelas and Deutch in the paper |^ . These authors consider 
a conformally coupled scalar field with Dirichlet or Neumann boundary conditions and electromagnetic field. Here we 
generalize the corresponding results for (i) an arbitrary number of spatial dimension, (ii) for a scalar field with general 
curvature coupling parameter and (iii) nonzero mass, (iv) for mixed boundary conditions of Robin type. In addition, 
(v) we consider the VEV's in the region between the plane and Rindler horizon for scalar and electromagnetic fields. 

To obtain the expectation values for the energy-momentum tensor we first construct the positive frequency Wight- 
man function (this function is also important in considerations of the response of a particle detector [Q). For the 
region on the right of the plate, ^ > a, the application of the generalized Abel-Plana formula to the mode sum over 
zeros of the function Ki^^ allows us to extract the purely Rindler part without boundaries and to present the additional 
boundary part in terms of strongly convergent integrals (formulas (p^), (p5|)). The expectation values for the EMT 
are obtained by applying a certain second order differential operator and taking the coincidence limit. First, in Sec. 
Ill, we consider the vacuum EMT in the Rindler wedge without boundaries. After the subtraction of the Minkowskian 
part the corresponding components are presented in the form ( ^l|) with notations (32). For the massless scalar the 
VEV's of the EMT take the form (|^). This expression corresponds to the absence from the vacuum of thermal 
distribution with standard temperature T — (27r^)~^. In general, this distribution has non-Planckian spectrum: the 
density of states factor is not proportional to uj^~^duj. The spectrum takes the Planckian form for conformally coupled 
scalars in D = 1, 2, 3. It is interesting to note that for even values of spatial dimension the distribution is Fermi-Dirac 
type. Further, on the base of procedure already used in |2^, for the massive case we present the vacuum EMT in 
another alternative form, formulas (p0|)-(p2|). By using the identity proved in Appendix]^ we show that for the cases 
of conformally and minimally coupled scalars these formulas are equivalent to those previously obtained in Ref. [p7| . 
For the massive scalar the energy spectrum is not strictly thermal and the corresponding quantities do not coincide 
with ones for the Minkowski thermal bath with temperature (27r^)^^. This fact is illustrated in Fig. |l|. In the limit 
3> 1 the vacuum EMT components are exponentially decreasing functions on m^. In the massless limit we obtain 
another equivalent representations, Eqs. ( p7|) and (|58|). 

In Sec. [lV| we investigate the boundary part of the vacuum EMT induced by a single plate in the region ^ > a, 
formula (|61|). The corresponding result for the massless ID field can be obtained by the analytic continuation and has 
form (|69|). In Appendix B| we show that this result coincides with the formula derived by direct calculation in D = 1. 
Further we investigate the asymptotic properties of the boundary EMT. Near the plate surface the total VEV's are 
dominated by the boundary part and the corresponding components diverge at the boundary. For non-conformally 
coupled scalars the leading terms are given by formulas ( [7l| ) and ( fz^ ) and are the same as for an infinite plane 
boundary in Minkowski spacetime with Dirichlet and Neumann boundary conditions. These terms do not depend on 
the mass or Robin coefficients, and have opposite signs for Dirichlet and Neumann cases. For large values of ^ a the 
leading term in the asymptotic expansion of the boundary induced EMT does not depend on the Robin coefficient B 
and has the form (^). In this limit (T;'=)(^V(r'/'>iub ~ ln"^(2^/a) and the vacuum EMT is dominated by the purely 
Rindler part. For a given ^ the boundary part of the EMT tends to zero as a — > 0, i.e. {Tj')'-''^ in the limit when 
the barrier coincides with the Rindler horizon. As a result the corresponding VEV's of the EMT are the same as for 
the FuUing-Rindler vacuum without boundaries. Hence, the barrier located at the Rindler horizon does not alter the 
vacuum EMT. 

The vacuum stresses in the region < ^ < a are investigated in Sec . for scalar and electromagnetic fields. The 
corresponding boundary parts can be presented as Eqs. (^Tj) and ( |107| ), respectively. These formulas differ from the 
ones for the region ^ > a by the replacements — > , K^^ I^^ . For the scalar field the boundary part diverges at 
the horizon with the leading behavior (|9^), (^). The divergence of the purely Rindler part is stronger and this part 
dominates near the Rindler horizon. Unlike to the scalar case, the boundary parts of the vacuum EMT components 
for the electromagnetic field vanish at the horizon (see asymptotic formulas (111)). In Fig. ^ we present the vacuum 
densities for the electromagnetic field, generated by a perfectly conducting plate in both RR and RL regions. As seen, 
unlike to the RR region, where the energy density is negative for all a < < oo, in the RL region the energy density 
is negative near the Rindler horizon, but is positive near the plate. As a result for some intermediate value of ^ we 
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have zero energy density. 
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APPENDIX A: SUMMATION FORMULA OVER ZEROS OF A'„ 

In Sec. H we have shown that the VEV's for the EMT components in the RR region contain sums over zeros 
Lo — LUn{r]), n = 1,2,... of the function Kit^jirf) = AKi^{ri) + bzKl^(r]) for a given 77. To obtain a summation formula 
over these zeros we will use the generalized Abel-Plana formula derived in |p^ , figt] . In this formula as functions f{z) 
and g{z) let us choose 

/(z) = — sinh7rzi^(z), (Al) 

TT 

9{z) = W~r\ 

with a function F{z) analytic in the right half-plane Rez > 0, and quantities with overbars are defined in accordance 
with Eq. (p^. For the sum and difference of these functions one has 

giz)±fiz)^2Fiz)^-p^. (A2) 

By using the asymptotic formulas for the Bessel modified functions for large values of index, the conditions for the 
generalized Abel-Plana formula can be written in terms of the function F{z) as follows 



\F{z)\ < ei\z\)e-'''' j , z^x + iy, x > 0, \z\ ^ 00, (A3) 

where |z|e(|z|) — > when \z\ — > 00. Let u — ujn be zeros for the function Kit^^rj) in the right half- plane. Here we will 
assume values of A and b for which all these zeros are real. It can be seen that they are simple. This directly follows 
from the relation 



KLA^V), (A4) 



which is a direct consequence of integral formula (|TJ) with — — Ki^^^{^ri). Substituting functions ( Al) 

into generalized Abel-Plana formula (formula (2.11) in Ref. [|l9| ) and taking into account that the points z — LUn are 
simple poles for the function g{z) and the relation liuj^l) — I-ii^niv): one obtains 



^™ ^ 1^ t Ma I 2 / smhTTzF{z)dz 



1 r iM 



27T Jo KM 



F(ze"/2) ^ F{ze-^''^)\ dz, (A5) 



where Un, < I < Um+i- Here we have assumed that the function F(z) is analytic in the right half-plane. However, 
this formula can be easily generalized for the functions having poles in this region. 
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APPENDIX B: L» = 1 MASSLESS CASE: DIRECT EVALUATION 



Here we show that the direct evaluation of the boundary VEV's in D ~ 1 gives the same resuh as the analytical 
continuation from the higher values of the space dimension. Now in equation (|l3| ) for the function <f>{£,) one has A = 0. 
The corresponding linearly independent solutions are e±*'^in(4/'>). The normalized eigenf unctions satisfying boundary 



spon 

condition (kW are in form 



e 



ifi^i^.T) = ^=ain[ujln{^/a) - P], < u < oo, (Bl) 

y'TTLO 

where 



A + iujB/a 



(B2) 



Note that these eigenfunctions have the same form for the regions ^ > a and ^ < a. Unlike to _D > 1 cases now the 
both types of t he e igenfunctions are bounded and the spectrum for oj in the region ^ > a is continous. Substituting 



eigenfunctions (Bl) into the mode sum formula after some algebra one finds the Wightman function in the form 
(P4|), where the first summand on the right is the D = 1 Wightman function for the Rindler wedge without boundaries: 

G+(x, x') = — cos [c. ln(e/e')], (B3) 

and for the boundary part one has 

1 /"OO 7 

{^{x)^{x'))^^^ =-— — e-'"("-"')cos[c^ln(ee7a')-2/3]. (B4) 
2tt Jo uj 

To transform this expression we write the cos function in terms of exponentials and rotate the integration contour 
by angle ±7r/2 for the term with exp[ja; ln(^^'/a^)] and by angle for the term with exp[— ja;ln(^^'/a^)]. Here 

and below the upper and lower signs correspond to the cases ^, ^' > a and 5, ^' < a respectively. Under the condition 
|r — t'I < I ln(5^'/c'^)l the integrals over the arcs with large radius in the complex w-plane tend to zero and we receive 

=-f r^4|^e--l-(«'/a^)leosh[c.(r-r')]. (B5) 
ZTT Jq a =f ujB/a 

The boundary contribution to the vacuum EMT can be found substituting this expression into formula (|^) . It can be 
easily seen that as a result for the boundary part we obtain formulas (^9|) and (|9^). Hence, we have shown that the 
direct evaluation gives the same result as the analytic continuation. 



APPENDIX C: PROOF OF IDENTITY (M) 



In this appendix we will prove identity (|5^), which was used to see the equivalency of our formulas (|50|)-(g2|) for 
the VEV's of the EMT in the special cases of the conformally and minimally coupled massive scalars to the results 
derived by Hill in Rcf. . Our starting point is the integral 

/•OO /"OO 

dkk^-^ dujij^e-'^^K^^XOK^^iXO (CI) 



with A defined as in Eq. (p^). By using formulas (^) and |24 

f (.v^) ^ 2<-=.'^,„-T (£^) (C2) 

for this integral one obtains 
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where 7 is defined by Eq. (ES). In particular taking = C we receive 



D-1 



dy 



(7r2+y2)3 ;,(D-l)/2 ' 



with z defined in accordance to Eq. (p9|) . Another form for can be obtained by using the relation 



d( 



The substitution into Eq. (CI) yields 

/:(e,^') = 2(^-3)/2^^D + lp 
1 - 



D- 1 
2 

1 d 



1^ ^ (^^)(D+l)/2 



d{'mS,)^ d{m^) 
For the coincidence case ^' = C from here one obtains 

'D~l\ r°° dy 



K(D-i)/2{m'-f) \ 
(m7)(^-i)/2 J 



„ „ , ,L> + 1 - (D - 1) coshy — ,„ iwo 



-(1 — cosh y) 



K(^D-l)/2{z) 



(C4) 



(C5) 



(C6) 



(C7) 



Comparing expressions (C4) and (|C7|) one obtains identity (B 
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